Trapped two-dimensional condensates with synthetic spin-orbit couphng 
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We study trapped 2D atomic Bose-Einstein condensates with spin-independent interactions in 
the presence of an isotropic spin-orbit coupling, showing that a rich physics results from the non- 
trivial interplay between spin-orbit coupling, confinement and inter-atomic interactions. For low 
interactions two types of half-vortex solutions with different winding occur, whereas strong-enough 
repulsive interactions result in a stripe-phase similar to that predicted for homogeneous conden- 
sates. Intermediate interaction regimes are characterized for large enough spin-orbit coupling by an 
hexagonally-symmetric phase with a triangular lattice of density minima similar to that observed 
in rapidly rotating condensates. 



Introduction. The engineering of synthetic electro- 
magnetism in ultra-cold gases has recently attracted a 
major attention [l|. Although an homogeneous effective 
magnetic field may be generated by simple rotation, re- 
cent techniques based on appropriate laser arrangements 
and tailored dressed states allow for more flexible con- 
trol of effective artificial gauge fields [3, Q • Interestingly, 
these techniques allow as well for the creation of non- 
Abelian gauge fields M JS, and more specifically spin- 
orbit couphng (SOC) [a-tlOj, a crucial effect in solid-state 
physics, essential for topological insulators [ll|. Recent 
ground-breaking experiments have explored this fascinat- 
ing possibility, demonstrating the creation of SOC with 
equal Rashba and Dresselhaus strengths (l^ . 

The creation of SOC in spinor gases opens fascinat- 
ing questions about the physics of ultra-cold gases with 
SOC, which have aroused a rapidly-growing theoretical 
attention both in what concerns degenerated fermions M, 
H, [l3| and Bose-Einstein condensates(BECs) 0, [l5l - l22f . 



In particular, Wang et al. have recently shown that the 
ground-state of an homogeneous, i.e. untrapped, two- 
component BEC with SOC is a single plane-wave phase 
or a spin stripe phase depending on spin-dependent inter- 
actions [l8| . On the other hand, non-interacting trapped 
spin-orbit coupled BECs are expected to present a half- 
quantum vortex configuration 0, ■ 

In this Letter, we show that the interplay between trap 
energy, SOC and interactions leads to a rich ground-state 
condensate physics. We consider in particular the com- 
plete phase diagram of trapped two-dimensional conden- 
sates with spin-independent interactions in the presence 
of isotropic spin-orbit coupling. This phase diagram is 
characterized by four different phases. As expected, for 
low interactions, a half-vortex solution (IIV(l/2) below), 
with angular momentum I such that \l + 1/2| = 1/2, is 
recovered. Increasing interaction leads to a second half- 
vortex solution (HV(3/2)) with a higher \l + 1/2| = 3/2. 
For a sufficiently large SOC the IIV(3/2) phase enters for 
larger interactions into an hexagonally-symmetric phase 
characterized by the appearance of a triangular lattice of 
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FIG. 1. Total density for «: = 15 and g = 0.05 (a), g = 
0.1 (b), 3 = 0.85 (c), and g = 2.0 (e). Figure (d) depicts the 
momentum distribution for g = 0.85, and Fig. (f) shows the 
spatial distribution of component 1 for g = 2.0. 

minima in the total density, similar to a vortex lattice in 
fastly rotating BECs ^23J . Finally, for even larger interac- 
tions, a spin-stripe phase develops, similar to that found 
in homogeneous BECs • We discuss the existence and 
main properties of these phases. 

Model. We consider a BEC of atoms with mass to, 
confined in two-dimensions on the xy plane, by a tight 
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harmonic confinement along z. The atoms have two avail- 
able internal states, which in typical experiments are, as 
mentioned above, atom-light dressed states of Zeeman 
sublevels with a proper laser configuration. Typically, 
inter-atomic interactions between Zeeman components 
are very similar. Hence, for simplicity we consider spin- 
independent interactions characterized by the 2D cou- 
pling constant g > 0, which hence determines as well the 
interactions between the atoms in the dressed states. 

Both components are equally confined in an isotropic 
harmonic potential V{t) = muj'^r'^ with r'^ = -\- . 
We consider that an appropriate laser configuration is 
chosen such that the atoms experience an isotropic gauge 
field A = hk{(Txex + (Jyey), where ax,y are the Pauli 
matrices and ex,y are the unit vectors along the directions 
X and y [f| . The physics of the condensate in the presence 
of SOC and interactions is described by the Hamiltonian 



2m 2 



*, (1) 



where Vl/ is a two component spinor. 

Homogeous solution. In the absence of trapping, 
V{r) = 0, the condensate minimizes the energy by spon- 
taneously breaking the rotational symmetry, develop- 
ing a__jj)atial spin modulation along an arbitrary direc- 
tion [18]. Chosing that direction as a;, we may re- write 
* = (^i{f^x)<^^(^^ For a constant * (with = n), both 
the interaction energy, gn, and {{~ihV — A)^) = 0, ac- 
quire its minimum value. Chosing = y/n{l,0), the 
spin modulation acquires the form: 



cos(kx) 
i sin kx 
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The presence of the trap modifies and enriches the 
ground-state physics since it introduces an additional 
energy scale, huj. The trapped BEC is hence best de- 
scribed by two dimensionless parameters, k = kIho, with 
= h/rrujj, and g = gm/ti^. Note that due to the 2D 
nature of the problem, the dimensionless coupling con- 
stant g is actually independent of uj. Below we employ 
dimensionless expressions, using oscillator units. 

HV(l/2) Phase. The non-interacting case, g = 0, 
is best described in momentum representation [7[. In 
the absence of trap, the Hamiltonian acquires the form 
(fik — A)^/2, which presents two eigen-energy branches, 
e_j_(fc) = (fc ± k)^/2, with k = |k|. The corresponding 
eigen-spinors are of the form 
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with if the polar angle of the vector k. Assuming a 
dominant SOC, such that k^/2 is larger than any other 
energy scale in the problem, we may reduce the anal- 
ysis to the lowest energy branch: ^(k) = ?/'(k)u_ ((ys). 
Introducing the trapping potential, V 



Vfc the gradient in momentum space, and using V'(k) = 
Y,i k~^/^fi{k)e'''^, one obtains: 

I / (? + 1/2)2 X /;^_~)2 

Note that the original kinetic energy term becomes 
a mexican-hat-like "potential" in momentum space, 
whereas the trap results in a " radial kinetic energy" term 
and a "centrifugal barrier" 0]. For k ^ 1, we may de- 
velop around k k, obtaining the eigenenergies 
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where ri characterizes the radial excitations of the 
mexican-hat. 

Without interactions the lowest energy is given by n = 
and |/ + f/2| = 1/2, which are states of the form: 
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These states present a spatial dependence of the form: 
*5'(r = (r,a)) ~ (Jo(Atr), e*" Ji(Kr)), and *^i(r) ^ 
(— e~^" Ji(ti:r), Jojur)), and hence constitute half-vortex 
solutions [I5|. In the following we denote this phase as 
HV(I/2). Solutions with / = and / = — I are degener- 
ated, being eventually selected by quantum fluctuations 
beyond our mean-field approach. However, any linear 
combination of / = and / = — I produces the same 
rotationally symmetric total density profile (Fig. [Ha)). 
The selection of Z = or / = — I would lead to a ring- 
like momentum distribution with fc ~ k, whereas linear 
combinations result in a cos ip modulation along the ring. 

HV(3/2) Phase. For ^ 1, angular excitations 
along the ring are much less energetic than radial ex- 
citations. Hence weak interactions result in the popu- 
lation of higher \l + 1/2 1 values. Simulations of Eq. ([T]) 
show that the HV(l/2) phase remains the ground state 
for sufficiently small g (Fig. [Ha)). However, the total 
density of the HV(I/2) solution has a maximum at the 
trap center due to the Jo{Kr) contribution. On the con- 
trary, the solutions with n — 0, \l + I/2| — 3/2, depend 
as J;=-2,i(Kr) and J;=-i,2(«^''), and hence have a mini- 
mum of the total density at the trap center, presenting a 
reduced interaction energy compared to HV(I/2). Solu- 
tions with I / -I- 1/2 1 = 3/2 become energetically favorable 
at g = gcr, when the interaction energy balances the en- 
ergy difference between both solutions for g = 0: 
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with /a = / 2rdre-°" 
2.35(277/^2). 



./i(l,At) /i(0,k) 
[j2„(Kr) + j2„+i(Kr)]". For large 
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-V^/2, with 



Decomposing ipCk) = evaluate ^ = 

(|/ + 1/2|) =^|az|2|Z + l/2| (Fig.[i;a)). For g < the 
phase HV(l/2) is characterized by a plateau at ^ = 1/2. 
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FIG. 2. (a) Value of ^ = {\l + l/2\) as a function of g for 
K = 15, the dashed lines indicate the value 1/2 and 3/2 char- 
acteristic of HV(l/2) and HV(3/2), respectively; (b) Aspect 
ratio X the cloud for the same case. The dashed curve 
indicates the expected aspect ratio from Eq. (|10p . 

At g — gcr a sudden jump occurs into ^ = 3/2. Inter- 
estingly, the numerical simulation of Eq. ([T|) shows that 

(2) 

up to a second critical gcr , a second plateau at ^ = 3/2 
occurs. We denote this phase as IIV(3/2), which is char- 
acterized by a ring-like density profile (Fig. [Ijb)). As for 
the HV(l/2), I = 1 and I = —2 are degenerated, and 
any linear combination of them produces the same ro- 
tationally symmetric total density profile. Again, as for 
HV(l/2), the selection of / = 1 or / = —2 leads to a 
ring-like momentum distribution whereas combinations 
of them lead to a cos(3(^) modulation along the ring. 

(2) 

Lattice phase. For g > gcr , C departs from the 3/2 
plateau, rotational symmetry is broken, and a differ- 
ent hexagonally-symmetric phase occurs, characterized 
by the developing of a triangular lattice of minima in 
the total density profile (Fig.[TJc)), which resembles that 
observed for vortex lattices in rapidly-rotating conden- 
sates [2^. The spacing between the minima is indepen- 
dent of g, being ~ 47r/3K, such that the border of the first 
Brillouin zone of the lattice lies on the k — k ring. On 
the contrary, as discussed below, the width of the cloud 
envelope a/ (x^) is independent of k for dominant SOC, 
being only dependent on g. Hence enhancing interactions 
just increases the number of observed density minima 
(oc K^(x^)) keeping invariant the lattice structure. The 
individual components present in typical numerical simu- 
lations an involved density and phase distribution, char- 
acterized by vortices and anti-vortices of different quan- 



tizations. In this Letter we are not interested in them, 
since the particular distribution among the components 
may depend on details of spin-dependent interactions. 

An interesting insight about the lattice phase may be 
gained from its momentum distribution characterized by 
the appearance of six maxima at angles ipj ~ i7r/3, 
around the ring of radius k (Fig. [I{d)). The lattice 
hence results from the combination of three pairs of op- 
posite momenta. Major features of the overall density 
profile may be obtained from these six-peaked structure 
by means of a Gaussian ansatz of the form: 



V'(k) = ^aje" 
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with aj = l/y/6, h 



kx cos ipj 



(8) 
and 



v<p = — fc^sin(pj 

„2\ 



ky cos fj . 



This solution is char- 
acterized by {x'^) = (y^) — {If + l'^)/4:. For a suf- 
ficiently large n, the energy may be approximated as 
2(x2) = y{l + ff){l + gf{r,)/3TT), with 
377'*, and 77 = lip/lr- Mini- 
•^^ in excellent agreement 



fir,) = I/77 + 8/V3 + IO7724 
mizing the energy we obtain [x~ 
with our numerical simulations (since the orientation of 
the cloud is arbitrary, the x direction is determined nu- 
merically as that with the larger width). Interestingly, 
the kinetic energy Ekin = 1/4/^ just contains the radial 
width Ir, since the dependence is exactly cancelled by 
the SOC. This cancellation leads to a global shrinking of 
the cloud for the lattice phase, and plays a major role in 
the stripe phase discussed below. 

Stripe phase. All the phases mentioned above ful- 
fill = {y^)/{^^) — 1; either due to a full rota- 
tional symmetry or due to hexagonal symmetry in the 
case of the lattice. However, at a sufficiently large g, 
the ratio x departs from 1 (Fig. [2jb)), i.e. the cloud 
breaks the hexagonal symmetry, becoming elongated 
along X (Fig. [BJe)). After a transient, in which x de- 
creases abruptly (Fig. [SJ^b)), the condensate acquires for 
a sufficiently large g a stripe form along the major axis 
X (Fig. mf)), similar to that obtained in the case of ho- 
mogeneous condensates with SOC [l8l |. 

The properties of the trapped stripe, and in particular 
its elongation along the direction of the stripe modula- 
tion, may be understood from the ansatz ([8]), but just 
considering two peaks on the ring at ipj = 0,7r. In this 
case, (x^) = and (y^) = 1^/2. For a sufficiently large 
K the energy functional may be approximated as 
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As mentioned above, the kinetic energy, which in absence 

reduces to 775-, since the 



of SOC is of the form -rri 

4lr 



1 



Icp contribution is exactly cancelled by the SOC. This 
breaks the symmetry between x and y axis, leading to 
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FIG. 3. Phase diagram as a function of k and g. The dashed 
line indicates the result of Eq. ((7|. 

an effective compression along y. Energy minimization 
leads to an equation for the aspect ratio 
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which provides an excellent agreement with our numeri- 
cal calculations as shown in Fig. [^Jbottom). 

Phase diagram. We have obtained the ground-state 
phase diagram by numerically solving Eq. ([1]) for a large 
number of g and k values. Our results are summarized in 
Fig- El The border g^cr between the lattice and the stripe 
phase is determined as the g value at which x — 0.75. 
Due to the I/k^ dependence of the angular modes, for 
large k values the HV(l/2) and HV(3/2) shrink, being 
confined to a regime of progressively smaller g values as 
K increases. On the contrary, the lattice phase becomes 
progressively more favorable extending for very large k 
almost to 5 = 0. The border between lattice and stripe 
approaches g — 1 for large k. When k decreases, the 
lattice phase shrinks until disappearing below Kc ~ 7.5. 
Below Kc there is a direct HV(3/2) to stripe transition, 
at which the rotational symmetry is broken. 

Conclusions. In summary trapped spin-orbit coupled 
BECs present intriguing ground-state properties as a re- 
sult of the interplay between trap energy, spin-orbit cou- 
pling and interactions. The corresponding phase diagram 
includes two rotationally symmetric phases, HV(l/2) and 
HV(3/2), a hexagonally-symmetric phase with a trian- 
gular lattice of minima, which resembles a vortex lattice, 
and a stripe phase, similar to that predicted for homoge- 
neous BECs. In experiments the phases should be clearly 
distinguisable either by monitoring in situ the density 
profile, or by time-of-flight imaging. The latter maps the 
momentum distribution, and, in particular, the lattice 
phase will be characterized by six out-going clouds. 

Finally, let us point that in our paper we have consid- 
ered for simplicity spin-independent interactions. In typ- 
ical experiments spin-dependent interactions are much 



weaker than spin-independent ones, and hence should not 
affect the overall density profiles, although the particular 
distribution among the components will change. How- 
ever, due to the dressed nature of the states in actual ex- 
periments, spin-dependent interactions do not simply re- 
duce to a different inter-component density-density inter- 
action, as in usual two-component condensates. On the 
contrary, phase-dependent terms will occur [l^], which 
may considerably complicate and enrich the physical pic- 
ture. This issue, and anisotropic effects, both in the trap 
and in the SOC, will be the focus of further investigation. 

Note added: When completing this paper we became 
aware of two recent works in which a similar problem 
is considered. In Ref. 2J|, which considers the effect 



of spin-dependent interactions, phases IIA (and lA), lA' 
(and part of IIB) and IB (and part of IIB) correspond 
to the IIV(l/2), Iiy(3/2) and lattice phases. A new re- 
cent version of Ref. [15| discusses the border of the stripe 
phase at small k < 4, showing a direct half- vortex to 
stripe transition, in agreement with our results. 

We acknowledge support from the Center for Quantum 
Engineering and Space-Time Research QUEST. 
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